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Moist stratified turbulence is studied in a two-dimensional Boussinesq system influenced 
I— I by condensation and evaporation. The problem is set in a periodic domain and employs sim- 

^ ^ pie evaporation and condensation schemes, wherein both the processes push parcels towards 

Q saturation. Numerical simulations demonstrate the emergence of a moist turbulent state 

^ consisting of ordered structures with a clear power-law type spectral scaling from initially 

^ spatially uncorrelated conditions. An asymptotic analysis in the limit of rapid condensation 

and strong stratification shows that, for initial conditions with enough water substance to 
^ ^ saturate the domain, the equations support a straightforward state of moist balance char- 

acterized by a hydrostatic, saturated, vertically sheared horizontal flow (VSHF). For such 
initial conditions, by means of long time numerical simulations, the emergence of moist 
balance is verifled. Speciflcally, starting from uncorrelated data, subsequent to the devel- 
opment of a moist turbulent state, the system experiences a rather abrupt transition to a 
regime which is close to saturation and dominated by a strong VSHF. On the other hand, 
initial conditions which do not have enough water substance to saturate the domain, do not 
T-H attain moist balance. Rather, the system remains in a turbulent state and oscillates about 

^ moist balance. Even though balance is not achieved with these general initial conditions, 

the time scale of oscillation about moist balance is much larger than the imposed time scale 
^ of condensation and evaporation, thus indicating a distinct dominant slow component in the 

moist stratifled two-dimensional turbulent system. 
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I. INTRODUCTION 

Water, in its different phases, plays a significant role in shaping climate of the Earth [1],[2]. From 
a dynamical perspective, the energy release associated with changes in phase of water substance has 
a profound effect on specific atmospheric phenomena as well as the large-scale general circulation of 
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the tropical and extra-tropical troposphere [3], [4]. To develop an understanding of the interaction of 
advection and condensation/evaporation, we study a moist stratified two-dimensional (2D) system 
which consists of nonlinear gravity wave-modes (including zero frequency vertically sheared hori- 
zontal flows - VSHF), and a dynamically consistent heat source due to condensation/evaporation 
(note that the kinematic advection condensation problem, where the advecting flow is prescribed is 
also of significant practical interest [5]). The decomposition of the stratified Boussinesq equations 
into vortical and wave modes is specified in Lelong and Riley [6] (applications can be found in 
Bartello [7], Sukhatme & Smith [8]), mathematical details and a description of the VSHF can be 
found in Embid & Majda [9] and Smith & Waleffe [10]. 

The stably stratified 2D Boussinesq equations (without the influence of water substance) have 
been studied as a model of thermal convection [11]. Numerical work on the scaling of kinetic 
and potential energy in this system is seen to be quite non-universal : in addition to the original 
Obhukov-Bogliano postulates, scaling laws for the temperature variance (or entropy) have been 
proposed and realized [12], [13], [15]. Unlike pure 2D turbulence, the stably stratified system does 
not conserve enstrophy. In spite of this simulations suggest the possibility of an inverse transfer 
of kinetic energy, the formation of a VSHF and a forward cascade of temperature variance (or 
potential energy) [13], [14], [15]. On the other hand, noting that the 2D system involves the nonlinear 
interaction of gravity waves, the susceptibility of these waves to parametric subharmonic instability 
leads to the anticipation of a robust transfer of energy to smaller scales [16], [17], [18]. Numerical 
work in a geophysical context has documented this forward transfer of energy [19], [20], [21], [22]. 
Further, the possible development of singularities [23] and fronts [24], [22] as well as the global 
well-posedness of these equations has also attracted significant attention [25], [26], [2 7], [28]. One of 
the main issues we investigate in the present work is whether the moist stratified two-dimensional 
system is capable of exploiting the available potential energy in a systematic manner so as to 
spontaneously organize spatially uncorrelated initial data. 

Two-dimensional stratified models have been used quite extensively in atmospheric settings, 
for example the tendency to organize condensing convection by gravity waves in 2D has been 
noted by Huang [29]. Similarly moist 2D models have been used for studying tropical cyclone 
formation [30], [31], squall-line thunderstorms [32], [33], atmospheric density currents [34], unstable 
modes in moist convection [35] as well as in elucidating properties of precipitating convection 
[36]. In another vein, 2D models arise in the so-called super-parameterization schemes for general 
circulation modeling wherein a high resolution 2D model is embedded in a coarser three-dimensional 
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(3D) model to help resolve features of geophysical systems which support an extended range of 
scales [37], [38], [39]. Quite naturally, as the 2D model feeds into the larger scale 3D flow, it is 
important to develop a feel for 2D moist turbulence in its own right. 

In addition to the scaling and structure formation, we hope to address the issue of balance and 
"slowness" in this simplified moist system. Indeed, the notion of balance and the existence (or not) 
of a slow-manifold in geophysical systems has a rich history [40], [41], [42]. We refer the reader to the 
text by Daley [43] for an introduction and to Vanneste [44] (and the references therein) for recent 
progress in the subject. Loosely, given the wide range of temporal scales inherent in geophysical 
systems, it is natural to inquire whether physically relevant information is actually contained in 
an independent simpler reduced dynamical system that evolves on a comparatively slower time 
scale. Usually such a simplified or reduced description is possible in certain limiting scenarios, 
the celebrated example being the emergence of quasigeostrophic dynamics in the limit of small 
Froude and Rossby numbers for a rotating and stratified fluid [9], [45]. In particular, for the system 
at hand, by imposing strong stratification (which implies the presence of high frequency gravity 
waves) and rapid condensation/evaporation (which implies rapid heating/cooling) we would like 
to see if it is possible to asymptotically extract an emergent state of balance or a reduced system 
which is able to capture the main features of the entire system. If in fact the analysis points to 
such a simplified description, it is also of interest to examine whether general initial conditions are 
capable of attaining this state. 

We begin by introducing the basic equations and the schemes employed for condensation and 
evaporation. The physical problem at hand is to consider the fate of perturbations to a basic state 
which is saturated and in hydrostatic balance. A linear analysis in a periodic domain reveals that 
the most unstable modes have a columnar structure whose growth rate has a particularly simple 
form in the limit of rapid condensation (and evaporation). The predicted emergence of columnar 
structures with appropriate growth rates is then realized in linear and non-dissipative simulations. 
Another aspect of the linear analysis is the presence of a distinct slow mode, this mode is also ob- 
tained via an asymptotic expansion of the suitably nondimensionalized nonlinear equations. The 
slow mode represents a state of moist balance and is characterized by a hydrostatic, saturated 
VSHF. To test the accessibility of this state of balance as well as scaling and structure formation 
in this moist system we then proceed to nonlinear simulations under conditions of strong stratifica- 
tion and rapid condensation/evaporation. Employing spatially uncorrelated initial conditions, we 
examine the development of a turbulent state which includes the formation of coherent structures 
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and the spectral scaling of the vapor, liquid, velocity and temperature fields. After the emergence 
of a turbulent state, initial conditions which have sufficient total water content to saturate the 
domain are seen to abruptly transition to the aforementioned state of moist balance. For initial 
conditions which are insufficient in the above mentioned sense, the flow exhibits similar scaling 
and structure formation. However, due to the lack of sufficient water substance, the system cannot 
achieve moist balance. Rather, our numerical examination shows that the system remains in a 
turbulent state and oscillates in an aperiodic or quasiperiodic manner about the state of moist 
balance. Finally, a summary and discussion conclude the paper. 



II. GOVERNING EQUATIONS 



The moist Boussinesq framework has been the subject of numerous recent investigations : for 
example, in the scaling of moist turbulence [46], the modelling of precipitation fronts [47], moist 
Rayleigh Bernard convection [48], studying droplet growth in a 2D setting [49], tropical cyclone 
formation [31] and squall- line thunderstorms [33]. The governing equations are : 



Du Vp' , q .^.^ 

-m= - + i^)T'k, 

Dt po To 

DT' C 

Dqv _ n Dqi _ 
Dt ~ Dt ~ 

In the above u = {u,w) is the 2D incompressible velocity field while T',p' are the potential 
temperature and pressure perturbations. C is the latent heat of condensation and Cp is the specific 
heat. We consider an equation of state for an ideal gas wherein p = pRT. In effect these equations 
isolate the effect of condensation (evaporation) via dynamical heating (cooling). The heating and 
cooling are represented by the source S which is described later in the text. We have not included 
a rain-out velocity or the influence of water substance on the buoyancy and the equation of state. 
A discussion of these effects and their role in the moist framework can be found in Klein & Majda 
[50]. 

As per the usual Boussinesq setup, the buoyancy frequency is given by {qX/Tq)^ where g, Tq, A 
denote the gravitational acceleration, reference potential temperature and the background potential 
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temperature gradient, respectively (i.e. T = Tq + Xz + T'). In addition, qv{x,z,t) and qi{x,z,t) 
are the mixing ratios of water vapor and liquid water, and qs{x,z) is a fixed saturation mixing 
ratio profile. In essence, the mixing ratios are governed by advection-condensation-evaporation 
equations, and q = Qv + Qi (the total water content) is materially conserved, i.e. Dq/Dt = 0. 

Physically, q^ is advected via the fluid flow and if it exceeds the saturation value at the new 
destination, the mixing ratio proceeds to relax back to its saturation value. Similarly qi is advected 
via the flow, and whenever qi > and q^ < qs, the parcel experiences evaporation. These statements 
are quantified via 



C = —{qy{x,t) - qs{x)) H{q^{x,t) - qs{x)), E = —qi H{qs{x) - q^{x,t)). (2) 

Tc Te 

E and C represent evaporation and condensation, and Tc, Te are the time scales associated with 
these two processes respectively. H{-) denotes the Heaviside step function. Equation (2) encodes an 
essential simplification that at any given instant of time a parcel experiences either condensation 
or evaporation. In the limit of rapid condensation (and evaporation), i.e. Tc,Te — )• 0, we have 
< qv{x,t) < qs{x) and < qi{x,t) < qi{xQ,0) + qs{xo) — where (xo,0) is the Lagrangian 
origin of the parcel at (x, t). One can see that the upper bound on q^ is local while that on qi is 
spatio-temporally non-local. Finally, the source in (I) is given hy S = C — E. 

Conservation properties : Apart from the material invariance of the total water substance 
{Dq/Dt = 0), (I) also preserves the material conservation of moist static energy given by -^(T + 
^qv) = 0. With regard to an energy conservation law, (I) with periodic boundaries yields 



|/[„^ + .^ + ^T'V^|/5r'. (3) 

Therefore, the heating and temperature fluctuations have to be correlated for a growth in per- 
turbation energy [3] — a detailed discussion of the conservation properties of moist systems that 
include bulk cloud microphysics can be found in Frierson, Paulius & Majda [47]. 

Vorticity- stream formulation : We introduce a streamfunction (■(/'), with u = —dzipjW = dxip 
and vorticity uj = —A^. Further we consider qg = qsiz) = qo — f3z, i.e. the saturation mixing ratio 
is a linear function of z. Defining q^' — qv — qs, (1) can be re- written as 
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~Di ~ 9x ' 

DT C 

-m = -^^^ + 

^ = /3V'.-5, ^ = 5. (4) 

Linear analysis : Linearizing (4) about our base state of rest, hydrostatic balance, no liquid 
water and saturation (i.e. S" = 0), we obtain a constant co-efficient system of partial differential 
equations. Substituting plane wave solutions of the form exp{[ifc -x — at]}, in addition to a distinct 
(7 = mode, we get 



where 61 = H{q^ — qg) and 82 = H{qs — q^). As expected, when {5i,52) = (1,0) i.e. during 
condensation, ^(3 > X leads to an unstable situation. Further, when Tc — )• (rapid condensation) 
(5) yields max(|cj|) = where 6 = — A and this maximal growth is attained for fields 

independent of z. Therefore, columnar or vertically coherent structures are the fastest growing 
modes in linear rapid condensation. 

The slow mode o" = is possible only when 5i = 62 = 0, i.e. qy = qs- In addition, this mode 
satisfies hydrostatic balance with w = implying u = u{z). Hence, the slow dynamics of the linear 
system consists of a hydrostatically balanced, saturated VSHF. 

Nondimensionalization : Nondimensionalizing the momentum equations and the incompress- 
ibility condition from (1), we obtain (dropping primes) 



Dw d<p Te 

du dw ^ . 

dx dz 

Here Eu = V /{poU'^) and Fr = U/{NL) are the Euler and Froude numbers, and (j) = p/po- L,V 
and U are the scale parameters of the length, pressure and velocity fields respectively. Tg is a nondi- 
mensional number defined by Tg = Q/{XL) where the potential temperature is nondimensionalized 
using the scale parameter Q. 
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In addition, the energy and water substance equations take the nondimensional form 
= -i^)^ + (7^)i — )[i'lv - qs)H{q^ - Qs) - qiH{q, - q^)], 

JJt i OpH Tc 

^ = i — )[-iQv - qs)H{qy - qs) + qiH{qs - qv)], 

^ = i — )[iQv - qs)H{qy - qs) - qiH{qs - (7) 

where Tq = U/L represents the advective time scale and qv,qi and qg have been scaled by Qq = 
max[(7s]. Therefore, taken together (6) and (7) are the nondimensional forms of the governing 
equations. 

Since our interest is in stratified systems, we set Fr = e < 1 in the following asymptotic 
analysis. Further, considering small potential temperature variations to the background state, we 
set Tg = e and choose Eu = 1. Clearly, the influence of heating on the dynamics is likely to be 
significant when the time scale associated with condensation (and evaporation) is smaller than the 
advective time scale, i.e. ^ = ^. This choice is also consistent with the time scales associated with 
advection (of the order of days) and condensation (of the order of a few hours) in the atmosphere 
[47] . Finally, given the large latent heat of water substance we scale (£Qo)/(Cpe) = \. With these 
substitutions, (6) and (7) read 



Du _ 84) 

~Di ~ ~~dx' 

Dw 86 1 ^ 

= + -T 

Dt 8z e ' 

w + -^[{qv - qs)Si - qiS2], 

8u 8w Q 
8x 8z ' 

^ = \[-{qv - qs)h + qi62], 

^ = ^[{qv - qs)h - qM- (8) 

Asymptotic analysis : As the linear analysis reveals an instability during condensation, we 
examine the dominant balances when 5i = 1 and 82 = 0. For small e, plugging in an asymptotic 
expansion of the form f = + e/^ + ... for all fields, at 0(l/e^) we obtain 



DT _ 1 1 
~Dt ~ 
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and at 0(l/e) we obtain 



, = -w'^ + ql. (10) 



Carrying forward to 0(1) we have 



o__d4f_ 



-w 



Dt dz 

dx dz ' 

0° 1 
Dt"' = 

Here D^/Dt = dt + u^dx + w^dz, i.e. it represents advection via the order zero velocity fields. 

Let us make a few observations from equations (9), (10) and (11). As = 0, the fourth 
equation of (11) gives = Qv- This establishes a hierarchical link between the nondimensional 
vertical velocity and vapor field, specifically = ql and = q'^. Further, adding the last two 
equations in (11) gives 



^(9S + 9?) = 0. (12) 

Therefore, the order zero water substance is materially conserved under evolution via the zero order 
flow. We now specialize to the case wherein qg = qsiz), keeping this in mind the vapor equation in 
(11) gives (using = ql) 



1 1 o^ 

If (13) is to hold for general saturation profiles, it implies that q}, = 0. But, at the present level of 
approximation, qv = q^^ ^qI- Therefore, = g<j and ql = can be true globally only if the initial 
water substance is sufficient to saturate the domain. 
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For such initial conditions, collecting the information in (9), (10) and the liquid water equation 
in (11), we see that the vapor field tends to saturation, and the excess vapor is converted to liquid 
water. Further, as = (when ql = 0), this implies a hydrostatic balance between (p^ and T^, 
and from incompressibility we obtain a VSHF, i.e. = u^{z,t). With this information, the liquid 
water evolution in (11) simplifies to 



As —7- Qs, we have 5i — )• and 82 = 0. Hence, this is the nonlinear version of the slow mode 
(cj = 0) encountered in the linear analysis of the problem. In fact, following Majda, Mohammadian 
and Xing [51], we refer to this as a state of moist balance. From the preceding analysis, we 
deduce that any initial condition with sufficient water substance may evolve, possibly at long 
times, towards moist balance. Whether such a state is actually attained is examined numerically 
in the latter portion of the present work. 

Quite clearly, initial conditions which do not have enough water substance to saturate the 
domain cannot satisfy = Qs and ql = globally, and therefore cannot achieve a state of moist 
balance. In fact, we intuitively expect such initial conditions to evolve around moist balance, with 
a continual interplay of condensation and evaporation in different portions of the domain. It is 
unclear if, under these circumstances, the gross features of the system will evolve at a timescale 
comparable to Tc (and Tg) or whether a slower timescale (for example, the advective scale) will 
come to dominate the dynamics. Indeed, another aim of our numerical work is to quantify the long 
term fate of such initial conditions. 



III. NUMERICAL EXAMINATION OF THE SYSTEM 



Numerically, we solve the system in the vorticity-stream form using a pseudo-spectral method 
in a periodic domain of size L^. The time stepping is via a fourth order Runge Kutta scheme. Dis- 
sipation is of the form of a hyperviscous damping given by (— l)"+^f„A"/ (for all fields /). In par- 
ticular, we adopt the formulation employed by Maltrud and Vallis [52] wherein f„ = /rms/ (^m^""^)) 
where km is the maximum resolved wavenumber. In all simulations the hyperviscous order is eight. 

The equations solved are (4), and we set g/To = X = f3 = 1/e (along with C/Cp = 1/e, this 
amounts to a solving (8) in vorticity stream form). The simulations are carried out for a range of e 
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values that spans IQ-^ to IQ-"^ for linear runs and 10 ^ to 10 ^ for the nonlinear runs. This gives 
the required leading order dependence of Fr ~ = e. The linear saturation profile is prescribed 
as Qs = Qo[l — z/L(i\, where Qo = L^/e. As /3 = A, the condition for linear instability reduces to 
C/Cp > 1. In the atmosphere, for water substance, C/Cp ^ 10^, therefore along with small Fr and 
Tg, the instability criterion is satisfied. In the nonlinear simulations we set this ratio to 1/e. For 
the linear runs, we have used smaller values of C/Cp to demonstrate the clear agreement between 
analytical and numerical growth rates. 

The simulations performed fall into three categories : (i) Linear simulations. These are non- 
dissipative runs performed with an intent to test the veracity of the linear analytical growth rate, 
(ii) Non-linear simulations. These fall into two groups (a) Short time high resolution runs to 
detect and quantify emergent coherent structures and their spectral scaling, (b) Long time lower 
resolutions runs to investigate the fate of differing sets of initial conditions, specifically those that 
have sufficient initial water substance and those that do not satisfy the sufficiency requirement. 



We begin with a set of linear non-dissipative simulations. Based on our rapid condensation 
(tc — )• 0) linear analysis, we expect vertically coherent structures to be the most unstable modes. 
As the maximal growth rate is independent of wavenumber, we examine the growth of spatially 
uncorrelated initial data with equal power at all scales. We choose Tc = Te = e = 0.1 and 
the resolution of the simulations is 512 x 512. Figure 1 shows the initial condition and its fate 
after linear evolution for 10 nondimensional units of time {t/r^. Clearly, in spite of the ambient 
stratification, in agreement with our linear analysis the emergent vapor field exhibits a degree of 
vertical coherence. 

The linear rapid condensation analysis, after substituting g/T^ = A = /3 = 1/e, yields a maximal 



growth rate of - y — 1. We proceed to compute the numerical growth rate of the total dry energy 
(LHS of equation 3) for varying e(= Tc = Tg) and C/Cp. In particular, we let e vary from 10~^ 
to 10~^ for C/Cp = 10,4 and 2 respectively. As seen in Figure 2, on appropriate rescaling of the 
time (i.e. plotting vs t/e), all curves for a fixed C/Cp collapse on top of each other. Further, the 
numerical growth agrees with the analytical estimate, and this agreement improves as e = Tc —)• 0, 
i.e. we approach the limit of rapid condensation. 



A. 



Linear evolution 
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FIG. 1: The first shows the scaled (i.e. q'y/Qo) spatially uncorrelated initial condition. The second panel 
shows the emergent scaled perturbed vapor field for Tc — t,, = e = 0.1 at t/e = 10 from the linear and 
non-dissipative simulation. 




FIG. 2: A plot of the log of the LHS of (3) Vs t/e for single scale initial conditions. Smaller e runs are 
carried out to progressively longer times, and for a fixed C/Cp the the curves collapse on top of each other. 
The dashed lines show the analytical scaling oi J — 1 for £/Cp = 10,4 and 2 respectively. 



B. Nonlinear evolution : moist turbulence 



We now proceed with nonlinear simulations to detect and quantify the emergence of scaling 
and coherence from initially spatially uncorrelated conditions. As before q/Tq = A = /3 = l/e, and 
^ = 1/e. The simulations are carried out at a resolution of 750 x 750 and we begin with e = 0.1. 

Figure 3 shows the evolution of the power spectrum of the vapor variance (i.e. (g(,)2). As the 
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initial condition is super saturated, at small times we have condensation over the entire domain 
which manifests itself in a drop in the total vapor variance. Subsequently, as observed via the 
sequentially numbered curves in the first panel of Figure 3, the spectrum evolves due a systematic 
growth of variance at all scales. As seen in the second panel of Figure 3, which shows the variance 
in selected wavenumbers as a function of time, after the initial period of decay due to condensation 
throughout the domain, the growth of vapor variance is almost exponential at all scales. Though, 
the growth begins earlier for progressively larger scales, hence in a short time the spectrum attains 
a characteristic negative slope. Finally, the variance at all scales levels off and the spectrum attains 
an invariant form (with fluctuations at the largest scales). 

To test the robustness of this behavior, we have repeated the simulations for smaller e, specifi- 
cally for e = 5 X 10^^ and 10~^. As mentioned, with decreasing e, we have stronger stratification 
and we also approach the limit of rapid condensation/evaporation. The invariant form of the vapor 
variance spectrum in these two cases is shown in the first panel of Figure 3. Clearly, the spectrum 
remains the same for the range of e considered. The growth of the vapor variance in different scales 
(not shown) also follows a similar pattern to that in the second panel of Figure 3. 

A similar growth of energy at all scales is also observed (not shown) in the velocity and tem- 
perature fields, though the invariant spectrum obtained follows different power laws. With regard 
to the kinetic energy we obtain a steeper spectrum (a power law with approximately -2 slope) 
while the temperature variance or potential energy exhibits shallower scaling (a power law with 
approximately -1 slope). Once again, these spectra show the same scaling for the range of e we 
have considered. 

Even though the various spectra show an evolution that is reminiscent of the inverse cascade in 
2D turbulence (though with a remarkable diversity of power laws for the various fields) , the mecha- 
nism of generation of this spectrum is quite different from the interscale transfer that characterizes 
2D flows. Specifically, in the present unstable setting, the vapor, velocity and temperature fields 
experience a growth of energy at all scales and the spectral form emerges due to the comparatively 
larger growth at large scales. 

Along with the spectral evolution of the various fields, we also notice the emergence of ordered 
structures from the initially spatially uncorrelated data. An example of this is can be seen in 
Figure 4 which, for e = 0.1, shows the perturbed vapor field in various stages of its evolution. We 
refer to this state, i.e. one with multiscale coherent structures, as a state of moist turbulence. In 
the following sections, we inquire into the long term fate of moist turbulence for differing classes 
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FIG. 3: First panel : A plot of the log of the vapor variance for e — 0.1 as it evolves in time. The numbers 
on the plots indicate the time ordering (note that the spacing between these time slices is not equal). The 
flat line is the initial power spectrum, which as prescribed is uniform in wavenumber. The uppermost 
two curves with dots and open circles show the invariant vapor variance spectrum for e = 5 x lO^'^ and 
e = 10^^ respectively. Second panel : The growth of variance in time with e = 0.1 for wavenumbers 
k — 10, 20, 40, 80, 120, 160 and 190, ordered from top to bottom. The vertical lines denote the times when 
the spectra are plotted in the first panel. 




FIG. 4: Evolution of a spatially uncorrelated initial condition for e = 10^^. The three panels show the 
emergent perturbed vapor field [q'^/Qo) for the time labels 4, 5 and 7 as referred to in Figure 3. The initial 
condition can be seen in the first panel of Figure 1. 

of initial conditions. 



C. Nonlinear evolution : long time behavior 

We now proceed to the long-time nonlinear evolution of the system. The simulations are carried 
out at a resolution of 350 x 350 for e = 10-\5 x lO'^ and lO'^. In all cases, the evolution was 
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seen to be similar. We focus on the case e = 10 ^, but comparative plots from other runs are also 
presented. 

1. Initial data with sufficient water substance 

We first consider initial data such that the water substance is sufficient to saturate the entire 
domain. To satisfy this condition we prescribe supersaturation everywhere, and as before (f^ is 
a uniformly distributed spatially uncorrelated field. As expected the vapor field evolves from 
its initially uncorrelated state, and coherent structures develop along with a systematic spectral 
evolution. In fact, the early behavior of the system and the development of moist turbulence is 
described in Section III B. 

After the spectra of the various fields attain an invariant form, we follow the evolution of the 
system. In particular, the first two panels of Figure 5 show the horizontal and vertical components 
of the kinetic energy as functions of nondimensional time for all e = 10~^ and 5 x 10^^. In both 
cases, at early times the components oscillate with comparable contributions to the kinetic energy 
from the vertical and horizontal portions. But, at a critical time [t^*) the system experiences a 
rather abrupt transition beyond which the horizontal component dominates the kinetic energy bud- 
get. Note that r^* decreases with e, i.e. the transition to horizontal fiow dominance is progressively 
more rapid as e — )• 0. Focussing on e = 10~^, Figure 6 shows the horizontal fiow after ^ = r^* ~ 5 : 
as is evident, we observe the formation of a VSHF. The VSHF in each case is broad and the fiow in 
each direction occupies approximately half the domain. Though not shown, the development of a 
VSHF is accompanied by a tendency to saturation. Further, as the vertical velocity becomes very 
weak, the system also experiences a progressively diminishing amount of heating from condensation 
and evaporation. In fact, as per (14), the liquid water and temperature are passively advected by 
the VSHF. 

Gathering these pieces together we argue that subsequent to the development of a moist tur- 
bulent state, at long times, with sufficient initial water substance, in accordance with the asymp- 
totic state of moist balance, the strongly stratified and rapidly condensing/evaporating system 
approaches a saturated VSHF. 
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t/E t/e 

FIG. 5: The first two panels show the contribution of the horizontal (blue) and vertical (red) components 
to the kinetic energy of the flow as a function of the nondimensional time (j) for e = 10^^ (first panel) and 
e = 5xl0~^ (second panel). 
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FIG. 6: The horizontal flow with sufficient initial water substance for e = 10^^ at | = 10. 
2. Initial data with insufficient water substance 

We now proceed to initial data where the initial water substance cannot saturate the domain. 
As these are very long time simulations, they are carried out at a resolution of 256 x 256 grid 
points. To satisfy this condition we supersaturate the domain in a localized Gaussian blob, while 
undersaturating it in a similar manner in the rest of the domain. A snapshot of the initial vapor 
field is shown in Figure 7. We ensure that q^{x,z,t = 0) > 0, the total initial water substance is 
less than the amount required for the domain to be saturated and that Qv' is a spatially periodic 
function. The region of supersaturation ensures that the system evolves without external forcing. 
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FIG. 7: The initial vapor field when the initial water substance cannot saturate the domain, the inset shows 
a cross-section of and Qs through the center of the domain. 

As with the previously considered sufficient class of initial data, the spectral scaling of the 
various fields evolves to the invariant state described in Section IIIB. After the development of 
a turbulent flow, we once again follow the evolution of the system. In particular, as in Figure 
5, the two panels of Figure 8 show the vertical and horizontal contribution to the kinetic energy 
as a function of nondimensional time for e = 10"'^ and e = 5 x 10^^. In contrast to initial data 
with sufficient water substance, clearly we do not observe a transition to horizontal flow dominance. 
Further, as expected, the domain does not tend to saturation and there is persistent heating/cooling 
throughout the domain. 

Examining Figure 8, we note that the large amplitude oscillations are much slower than the 
imposed scale of condensation/evaporation (this is evident as the time in the plots has been nondi- 
mensionalized by Tc). Figure 9 shows the variance preserving temporal spectrum of the horizontal 
and vertical components of the kinetic energy for e = 5 x 10^^. As expected, there is a large amount 
of power at small frequencies or large time scales. Further, we can clearly see multiple peaks that 
lie between Tc and Tc/e which indicates a slow (as compared to Tc) aperiodicity in the time series. 
Also, at small time scales (or large frequencies), especially smaller than Tc, the flattening of the 
red-compensated spectrum shown in the inset of Figure 9 indicates a red noise scaling. 

The broad nature of the peak in the inset of Figure 9, as well as the multiple peaks seen in 
the variance preserving form of the spectrum, suggests an aperiodic or quasiperiodic character to 
the slow oscillations seen in the components of the kinetic energy. Performing a spectral analysis 
of the domain averaged heating as well as the mean vapor and liquid water content (not shown) 
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FIG. 8: Same as Figure 5 but with insufScient initial water substance as shown in the second panel of 
Figure 6. The simulations have been carried out for a significantly longer duration than shown in these 
plots. 

reveals a similar broad peak at long time scales compared to Tc. Quantitatively, the slow time 
scale has a range from double to approximately 1/e times the imposed condensation timescale. We 
have verified that simulations with other values of e (specifically, e = 10~^ and 10"^) yield similar 
results. 

Therefore, as expected from the asymptotic analysis, when the initial conditions do not have 
enough water substance to saturate the domain, the system does not attain a state of moist balance. 
In fact, the system oscillates about moist balance. This can be understood from Figure 8 where we 
see that, as the contribution of the horizontal flow kinetic energy increases, the system approaches 
moist balance. But, as the entire domain can never be saturated, the vertical flow kinetic energy 
soon gains prominence. This cycle repeats aperiodically at the aforementioned slow time scale. 
It is important to note that while the domain averaged characteristics are modulated at a slow 
time scale, the various fields maintain their spectral scaling and the system remains in a state of 
moist turbulence. Finally, even though it is not possible to associate a single timescale with the 
slow oscillation, we see that its broad peak extends from twice to 1/e times the imposed scale of 
condensation (and evaporation). 

IV. CONCLUSION AND DISCUSSION 

We have considered a 2D moist Boussinesq system in a periodic setting with simple condensation 
and evaporation schemes. In particular, these schemes are designed to push individual parcels 
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FIG. 9: The power spectrum of the kinetic energy component time series seen in Fig. (8) for e = 5 x 10~^. 
We consider an average over twelve long segments, each encompassing multiple oscillations seen in the signal. 
The main plot shows the variance preserving form of the spectrum. The inset shows the red compensated 
spectrum, i.e. spectrum after a multiplication by the square of the frequency. 

towards a saturated state. This allows possibly the simplest scenario wherein we can study the 
interaction of advection with evaporation and condensation in a dynamically consistent manner. 
Linearizing the problem about a state of saturation, with no liquid water and hydrostatic balance 
showed that the most unstable modes were vertically coherent. In fact, in the limit of rapid 
condensation, the form of the growth rate of the unstable modes took a particularly simple form. 
A suite of non-dissipative linear simulations verified these predictions. 

One of the main questions we wished to address in this work was the development of initially 
spatially uncorrelated vapor fields. Keeping in mind that the 2D Boussinesq system does not 
conserve enstrophy, and that there is significant theoretical and numerical evidence showing that 
the system supports a transfer of energy from large to small scales [16], [17], [18], [19], [20], [21], [22], 
it was not a priori evident to us that initially uncorrelated fields would develop into coherent 
structures, or that we would see a growth of energy at large scales. Pursuing nonlinear simulations 
we observed that coherent structures do emerge spontaneously from spatially uncorrelated initial 
conditions. Further, the scaling of the vapor variance as well the kinetic and potential energy 
attained an invariant form and follows a power-law. In fact, the system attains a state that can be 
aptly described as moist turbulence. 

Interestingly, during the nonlinear simulations, the growth of energy at large scales cannot be 
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attributed to an inverse transfer kind of mechanism seen in pure 2D turbulence wherein energy 
is transferred from small to large scales. In fact, in an unstable setting, the variance in all fields 
grew very rapidly at all scales and this rapid growth persisted for a longer time at larger scales 
leading to a characteristic negatively sloped spectrum. Given that the emergence of coherent 

structures in the vapor field did not depend on the conservation of cnstrophy, we are optimistic 
that a similar development of multiscale structures from initially uncorrelated fields will be possible 
in a three-dimensional moist stratified system. 

Our other goal was to examine the notion of balance and slowness in physically relevant lim- 
its of this idealized moist system. Imposing strong stratification (implying the presence of high 
frequency gravity waves) along with rapid condensation (implying rapid heating/cooling via con- 
densation/evaporation), wc found that, if there is sufficient initial water substance to saturate 
the entire domain, then the system supports a hydrostatic and saturated VSHF (a state of moist 
balance). If the initial water substance is not sufficient to saturate the domain, then the state of 
moist balance is inaccessible and the system evolves via a continual interaction of evaporation and 
condensation in different regions of the domain. 

Wc began our numerical investigation of the long time development of the system by considering 
initial conditions that have sufficient water substance. Simulations were carried for for a suite of 
e values (for progressively stronger stratification and smaller condensation timescales), and the 
evolution in each case was similar. As per expectations the system evolved with coherent structures 
and scaling emerging from initially uncorrelated conditions. After the development of invariant 
spectra, we followed the horizontal and vertical components of the kinetic energy as a function 
of nondimensional time (t/e). At small times, we noted oscillations with equal contribution from 
either component at different times. However, at ^ = r^*, we observed a rather abrupt transition 
beyond which the horizontal component dominated the kinetic energy budget. Examining the state 
of the flow after | = t^* , wc noted the formation of a persistent and dominant VSHF. In addition, 
the vapor field tended to saturation for | S> Tg*. In essence, with sufficient initial water substance, 
we saw that, the system (at long times) tends to the asymptotic state of moist balance. 

We then proceeded to initial conditions which were insufficient in initial water substance. The 
now generic behavior involving coherent structures and scaling was observed, but as expected 
from the asymptotic analysis, the system did not evolve to a state of moist balance. In fact, 
an analysis of the vertical and horizontal kinetic energy (as well as other domain averaged time 
series such as the heating and potential energy), revealed dominant slow aperiodic or quasiperiodic 
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oscillations and a fine temporal structure within each such oscillation. The period associated with 
these slow oscillations was much larger than the imposed condensation/evaporation time scale (tc), 
and therefore the system was modulated on a relatively slow time scale. In particular, the slow 
oscillation was seen to range from twice to 1/e times the imposed timescale of condensation (and 

evaporation). We emphasize that the spectral scaling of the various fields (noted earlier) persists 
in these conditions, i.e. even though the system experiences a series of slow oscillations about the 
state of moist balance, all fields maintain their characteristic scaling and, in essence, the system 
remains in a state of moist turbulence. 
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